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USSON  Group  Programs 

1.1  Introduction 

The  POISSON  group  progress  are  a  set  of  computer  codes  that  solve 
both  electrostatic  and  magnetostatic  problems  using  a  deformable 
triangular  mesh.  This  set  of  codes  has  been  modified  to  solve 
the  space-charge  flow  problem. 

The  first  step  in  solving  a  problem  using  the  POISSON  codes  is  to 
construct  a  non-uniform  triangulated  mesh  which  represents  a 
discretisation  of  the  problem.  The  program  which  generates  the 
mesh  is  called  LATTICE.  LATTICE  naps  equilateral  triangles  in 
logical  space  into  a  non-uniform  triangular  mesh  in  real  space. 

Logical  space  consists  of  an  array  of  regular  equilateral 
triangles  composed  of  three  sets  of  straight  lines  intersecting 
each  other  at  60  degrees  as  shown  in  Fig.  1.1.  Bach  vertex  of 
a  triangle  has  sssociated  with  it  a  logical  coordinate  (K,L) 
where  E  is  the  column  order  and  L  is  the  row  order.  LATTICE  will 
assign  each  (K,L)  coordinate  within  the  boundary  to  an  (x,y) 
coordinate  in  real  space.  To  illustrate  this  mapping  from 
logical  space  to  real  space,  the  logical  diagram  for  two 
concentric  circles  is  shown  in  Fig.  1.2a.  For  each  logical 
coordinate  on  the  boundary,  a  physical  coordinate,  (x,y)  must  be 
assigned  to  it.  For  the  concentric  circle  problem,  the  logical 
coordinates  on  the  boundary  are  assigned  the  physical  coordinates 
of  a  circle.  The  process  of  mesh  generation  results  in 


distorting  triangle  sides  in  the  logical  mesh  to  fit  all  region 


j)  Logic  dlogrom  for  two  concentric  circles. 
>)  Relaxed  triangulated  mesh. 
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boundaries  and  determining  the  location  of  all  interior  nesh 
points.  The  final  result  of  the  nesh  generation  is  shown  in  Fig. 
1.2b.  Note  that  the  anount  of  triangles  is  conserved  fron  the 
transformation  from  logical  space  to  real  space,  but  they  are  no 
longer  equilateral. 


1.2  POISSON  Group  Programs  Used  to  Solve  the  Electron  Flow  Problem 

The  following  is  a  list  of  programs  with  descriptions  of 
their  function  in  solving  the  space-charge  flow  problem: 

1.  AUTOMESH  -  prepares  the  input  data  for  LATTICE  fron 
physical  geometry  data  describing  the  problem.  By 
defining  the  physical  region  boundaries  and  specifying 
the  nesh  size,  AUTOMESH  will  construct  the  "logical" 
path  in  the  mesh  and  generate  (x,y)  coordinates  on  the 
boundary  for  straight  lines,  arcs  of  circles  and 
segments  of  hyperbolas. 

2.  LATTICE  -  generates  an  irregular  triangular  mesh  from 
input  data  which  includes  "logical"  and  physical 
coordinates  describing  the  boundary  and  boundary 
conditions. 

3.  POISSON  -  solves,  by  successive  point  over-relaxation, 
Poisson's  equation  for  the  scalar  potential  two- 
dimensional  problems,  calculates  the  derivatives  of  the 
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potential  (fields  and  gradients),  calculates  the  stored 
energy  and  performs  haraonic  (aultipole)  analysis  of 
the  potential. 


1.3  Units  Used  In  AUTOMESH.  LATTICE,  and  POISSON 


All  calculations  in  the  POISSON  group  programs  are  done  in  MKS 
units  only.  The  user  may  use  units  other  than  MKS  for  length, 
but  aust  appropriately  set  the  scale  factor  to  the  correct  value. 
This  will  be  discussed  later. 


The  following  is  a  list  of  units  used  in  the  program: 


length 


current 


(meters)  or  user  defined  by  specifying  Con(9) 
not  equal  to  1.0. 

(ampere) 


scalar  potential  V  =  (volt) 


electric  field  B 


induction  B 


(volt/m) 


( tesla) 


1.4  Format  Free  Input  Routine 


The  POISSON  group  codes  have  a  format-free  input  subroutine, 
FREE,  which  appears  in  the  computer  code  as  follows: 


Subroutine  FREE  ( N , RAY1 , N1 , RAY2 , N2 , RAY3 , N3 ) 


where  RAY1,  RAY2,  and  RAY3  are  the  array  (and/or  variable)  names 
which  are  to  be  read  in  and  Nl,  N2,  and  N3  are  the  numbers  of 
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values  RAY1,  RAY2  and  RAY3  are  to  have  respectively.  The 
variable  N  indicates  how  aay  naaes  are  to  be  read  in,  i.e.,  for 
N=l,  then  only  RAYl(Nl)  will  be  read  in. 


When  the  program  requires  specific  data  to  be  read  in,  or  certain 
variables  to  be  nodified,  the  following  command  will  be  executed: 


CALL  FREE  ( N, RAY1 , Nl , RAY2 , N2 , RAY3 , N3 ) 


The  user  will  then  input  the  data  using  a  specific  format. 


Special  input  characters  and  their  functions  are  described  below: 


"+”  and  indicate  the  sign  of  the  number  (mantissa  or 

exponent).  The  "  +  "  sign  is  optional. 

The  decimal  point  is  required  for  all 
floating  point  numbers. 

"E"  The  signed  integer  number  following  the 

letter  "E"  gives  the  power  of  ten  to  which 
the  previous  floating  point  number  is  raised. 
No  blank  before  or  after  the  "E"  is  allowed. 

"*"I  Store  the  number  following  "I"  into  the  array 

location  (I).  Successive  input  numbers  will 
be  stored  into  array  locations  (I),  (1+1), 

(1+2),  etc. 

"R"N  Store  the  previous  input  number  "N"  times 

into  the  array  locations  (I),  (1+1),  (1+2),  . 
.  .  through  (I+N-l).  A  blank  separating  the 

previous  input  number  and  the  "R”  is 
optional . 

"S"  Skip  the  rest  of  the  Nl,  N2,  or  N3  values 

requested  in  the  call  and  go  on  to  the  next 
array  or  return  if  the  current  array  is  the 
last  array  in  the  argument  list. 


”C"  Skip  the  rest  of  the  Nl,  N2  or  N3  values 


-  6  - 


requested  in  the  call,  set  Nl,  N2  or  N3  equal 
to  the  nuaber  of  values  input  into  the 
current  array  and  go  on  to  the  next  array  or 
return.  When  using  this  feature,  be  sure 
that  the  Nl,  N2  or  N3  in  the  calling 
stateaent  is  a  variable  and  not  a  constant. 

"  "  and  Blank  and  coaaa  are  the  only  non-nuaeric 

characters  allowed  in  the  input  field,  and 
one  of  these  characters  Bust  be  used  to 
separate  input  values. 


Coaaents  aay  follow  the  last  HSn,  "C"  or  required  nuaber  on  any 
input  line  (card). 


An  exaaple  is  given  below  which  illustrates  aost  of  the  features 
in  the  foraat-free  routine.  The  following  is  a  calling  sequence 
of  FORTRAN  stateaents;  A(I)  and  B(I)  are  diaensioned  arrays  and  K 
is  a  single  variable. 


DIMENSION  A(5) ,  B(100) 

N  =  100 

CALL  FREE  (3 , A , 5 , B , N, K, 1 ) 

When  the  call  to  FREE  is  involved,  the  prograa  requires  input 
data  given  either  froa  the  terainal  or  input  file.  The  following 
line  is  typed  on  the  terainal  for  the  above  call  to  FREE: 


-3. 0,4.0  +5 . 3B-2  R2  S  *20  0.1R10  COMMENTS  HERE 


Once  this  line  has  been  read  in  A,B  and  K  take  on  the  following 
values , 


A ( 1 ) =-3 . 0 
A ( 2 ) =4 . 0 
A ( 3 ) =0 . 0053 


A(4)=0. 0053 

B(20)  through  B(29)=0.1 
K=13 

and  tha  returned  value  of  N  is  10. 

1.5  Program  Flow  Chart  of  the  POISSON  Group  Progress 

Fig.  1.3  shows  the  coaaunicat ion  between  the  POISSON  group 
prograas  listed  in  Sec.  1.2. 

The  boundary  generator,  AUTOMESH,  requires  the  data  file 
AUTOMESH . DAT  which  is  created  by  the  user.  AUTOMESH  will  create 
the  data  file  LATTICE.DAT,  which  is  read  in  by  LATTICE. 

The  aesh  generator,  LATTICE,  requires  the  data  file  LATTICE.DAT 
which  ia  either  created  by  AUTOMESH  or  the  user.  LATTICE  creates 
two  files: 

1.  LATTICE. OUT,  contains  information  on  the  aesh  generation. 

2.  POISSON.DAT,  contains  the  aesh  and  boundary 
inforaation  that  is  read  by  POISSON.  LATTICE  writes 
this  data  file  in  unforaatted  output,  therefore  it 
cannot  be  viewed  since  it  is  written  in  internal 
representation  (binary). 

The  Poisson  solver,  POISSON,  requires  the  data  file  POISSON.DAT. 
POISSON  creates  the  following  two  files: 

1.  POISSON. OUT,  which  is  written  on  unit  90,  contains 
information  on  the  solution  of  the  electron-flow 


problem . 


AUTOMESH.DAT 


GENERATE 
LOGICAL  - 
BOUNDARY 


AUTOMESH 


LATT1CE.DAT 


GENERATE  THE 
MESH 


_ j 

LATTICE.  OUT 


DUMP  0 


SOLVE  THE 
ELECTRON  FLOW 
PROBLEM 


POISSON.DAT 


— ■*  poisson  - y 

POISSON. OUT 


DUMP  >  0 


POISSON.DAT 


PLOT 


FI9.  1.3  Program  Flow  Chart  of  the  POISSON  Group  Programs 


2.  POSITION. OUT,  which  la  written  on  unit  91,  contains  the 

trajectory  inforaation  and  ia  optional. 

Tha  plotting  prograa,  PLOT,  raquiraa  the  data  file  POISSON.DAT 
which  is  read  fron  unit  95. 


The  prograa  AUTOHBSH  finds  the  path  in  the  logical  eesh  which 
best  fits  the  physical  boundary.  Once  AUTOMBSH  has  found  this 
logical  path,  the  physical  coordinates  (x,y)  are  assigned  to  each 
logical  point  (K, L) .  This  prograa  only  determines  the  logical 
coordinates  for  the  physical  coordinates  on  the  boundary.  The 
prograa,  LATTICE,  deteraines  the  physical  location  of  the  logical 
coordinates  within  the  boundary.  AUTOMBSH  can  be  thought  of  as  a 
preprocessor  for  LATTICE. 

The  algorithm  in  AUTOMBSH  which  constructs  the  logical  aesh  from 
data  describing  the  boundary  of  the  problea  can  be  visualized  as 
follows: 

1.  Plot  to  scale  the  physical  boundary  of  the  problea. 

2.  Choose  a  triangle  size  for  the  aesh. 

3.  Overlay  the  array  of  identical  triangles  onto  the 
physical  boundary  and  determine  the  logical  path. 

This  procedure  of  constructing  the  logical  aesh  is  shown  in  Fig. 


The  input  date  for  AUTOMBSH  ia  read  by  a  "pseudo”  naaelist 
routine.  There  are  three  groups  of  input  data  for  AUTOMBSH  which 
ia  deacribed  below  in  aectiona  2.3,  2.4  and  2.5.  Data  groups  2.4 
and  2.5  are  repeated  for  each  region.  The  values  of  the 
variables  for  the  firat  region  are  preserved  for  succeeding 
regions  until  they  are  changed  in  the  next  region. 

2.2  Foraat  of  the  Input 

In  the  original  AUTOMBSH  prograa,  the  input  data  was  read  in  by 
the  FORTRAN  "NAMBLIST"  routine.  Since  the  HP-1000  does  not 
support  NAMELIST,  the  progran  was  rewritten  with  a  "pseudo" 
naaelist  routine.  This  aodified  code  tries  to  enulate  the 
NAMELIST  feature  as  closely  as  possible.  The  liat  of  variables 
in  the  NAMELIST  specifies  those  variables  which  nay  or  nay  not 
have  their  values  defined  in  the  RBAD  atateaent.  Not  all 
variables  of  the  NAMBLIST  need  be  read  in  nor  nust  th^ey  be  in  any 
order . 

The  very  firat  iten  to  be  input  is  an  integer  constant  which 
represents  the  nuaber  of  input  variables  to  be  read  in.  The  next 
iten  or  iteas  are  the  input  variables.  Each  input  variables  is 
represented  by  a  pair  of  FORTRAN  variables.  The  first  FORTRAN 
variable  is  a  character  variable  which  represents  the  naae  of  the 
input  variable  and  the  second  variable  is  the  nuaerical  value  of 


the  input  variable.  Since  an  unforaatted  READ  ststeaent  is  used 


to  read  io  the  pair  of  variables,  the  character  variable  euat  be 
included  in  single  quotes  and  either  a  blank  or  a  coaaa  aust 
exist  between  each  variable. 

As  an  exaaple,  say  the  following  three  input  variables  are  to  be 
read  in,  or  changed  froa  their  default  valuea:  ’NREG=’,  ,XMIN=', 
*YMIN=’.  The  following  two  lines  would  be  typed  in  the  data  file 
AUTOMESH . DAT : 

3 

*  NREG- *  4,  *  XMI N- *  2.0,  * YMIN= '  1.0 

The  integer  constant  3  in  line  one  indicates  that  there  are  three 
input  variables  to  be  read  in.  The  FORTRAN  variable  "NREG"  will 
be  assigned  a  value  of  4  and  the  FORTRAN  variables  "XMIN"  and 
"YMIN”  will  be  equal  to  2.0  and  1.0  respectively.  Note  the  use 
of  blanks  and  coaaas  in  line  two.  This  is  required  since  each 
itea  is  a  piece  of  data  and  aust  be  separated  by  a  blank  or  coaaa 
when  using  a  unforaatted  input  stateaent.  The  data  iteas  that 
appear  in  the  input  file  need  not  appear  in  any  special  order  and 
blanks  aay  be  used  freely  between  data  iteas  to  aake  the  file 
easier  to  read.  Also  note  that  the  three  pair  of  variables  in 
line  two  do  not  all  have  to  be  on  the  sane  line  since  the 
execution  of  an  unforaatted  RBAD  stateaent  is  not  finished  until 


all  the  listed  variables  are  assigned  values. 
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The  first  data  line  can  have  anything  in  coluans  2  through  80. 


For  an  electron  flow  problea,  coluan  one  Bust  be  blank.  If 
coluan  one  is  non-blank,  then  this  data  file  is  for  a  SUPERFISH 
problea.  If  coluan  one  is  blank,  then  this  data  file  is  for  a 
POISSON  problea.  SUPBRFISH  is  a  separate  prograa  which  solves 
for  resonant  nodes  in  cavities.  The  characters  in  coluans  2 
through  80  are  stored  for  run  identification  in  the  printed 
output . 

2.4  Input  Variables  That  Describe  The  Mesh  and  Boundary 
Conditions 

The  following  is  a  reference  list  of  input  variables  that 
describe  the  physical  coordinates  of  the  boundary: 


NRRG, 

IRBG , 

NAT,  CUR,  DBN,  ITRI, 

IBOUND ,  DX,  DY 

,  XMIN, 

XMAX, 

YNIN, 

YNIN,  YNAX,  NPOINT , 

IPRINT,  XREG1, 

XREG2, 

YRBG1 , 

YRBG2 

,  KRBGl ,  XRBG2,  XMAX, 

LRBG1 ,  LRBG2 , 

LMAX 

A  coaplete  description  of  each  variable  is  given  below.  Default 
values  are  ahown  in  parentheses.  If  the  user  desires  to  use  a 
default  value,  then  onit  the  data  iteas  froa  the  input  file. 


*  NRIG  =  ’ 


The  nuaber  of  regions  to  be  read  in  this  run 
(input  only  once  in  the  date  file). 


*  1810=’ 


The  LATTICR  region  nuaber.  The  default  value  is 
one  for  the  1st  region  data  group  and  is 
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incremented  by  one  for  succeeding  region  data 
groups. 


*MAT=*  Material  code  for  the  region. 

(1) 

=0,  all  points  inside  this  region  are  oaitted. 
=1,  air  (mu=l) 

=  2,  iron  with  internal  preaeability  table 

=3 ,  iron  with  input  permeability  table  no.  1 

=4,  iron  with  input  permeability  table  no.  2 

=5,  iron  with  input  permeability  table  no.  3 


*CUR= ' 
(0.0) 


The  fixed  potential  value  (volts)  on  the  boundary 
of  the  region.  For  magnetic  problems,  it  is  the 
total  current  (amps)  in  the  region. 


*D1N=*  The  charge  density  (colouabs/length**2)  in  the 

(0.0)  region. 


*ITRI=’  Triangle  type  generated  in  the  region 

(0) 

=0,  equal  weight  (default) 

=1,  equilateral 
= 2 ,  right 


* IBOUND=’  Region  boundary  indicator  (See  LATTICE  section  for 

a  more  complete  description) 

=  -1,  fixed  potential  Dirichlet  boundary  condition 
=  0,  Dirichlet  condition  on  the  boundary  (magnetic) 

=  1,  Neumann  condition  on  the  boundary 


*DX=’  The  width  of  the  triangles  in  the  logical  mesh. 


*DY='  The  height  of  the  triangle  in  the  logical  mesh. 

Default  values  for  ITRI=0,  DY=DX*SQRT( 3 . 0) /2 . 0 ,  or 
DY  =  DX  if  ITRI=2. 
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• XMIN* * 

(0.0) 


' XMAX= * 


*  YMIN*  * 

(0.0) 


* YMAX=  * 


*  XRBG1=  ’ 


*  XREG2= ' 


* YREGl*  * 


*  YRBG2= * 


'  KREG1= ' 

*  KRBQ2= ’ 

*  KHAX= * 

' LRBG1= * 
*LREG2=’ 
* LMAX= ' 


* NPOINT= ’ 


Miniaua  X  value  of  the  problea. 


Maxiaua  X  value  of  the  problea. 


Miniaua  Y  value  of  the  problea. 


Maxiaua  Y  value  of  the  problea. 


The  width  of  the  triangles  in  the  logical  aesh 
will  be  set  to  approxiaately  2.0  *  DX1  for  X 
values  greater  than  XRBGl. 


The  width  of  the  triangles  in  the  logical  aesh 
will  be  set  to  approxiaately  2.0  *  DX1  for  X 
values  greater  than  XREG2. 


The  height  of  the  triangles  in  the  logical  aesh 
will  be  set  to  approxiaately  2.0  *  DY  =  DY1  for  Y 
values  greater  than  YRBG1. 


The  height  of  the  triangles  in  the  logical  aesh 
will  be  set  to  approxiaately  2.0  *  DY1  for  Y 
values  greater  than  YREG2. 


When  the  ’’doubling”  scheaes  set  up  by  the  XREGn 
and  YREGn  variables  are  not  appropriate,  regions 
in  the  logical  aesh  with  different  triangle  sizes 
aay  be  explicitly  defined  by  inputing  the 
variables  KREGn  and  LREGn.  With  these  variables, 
one  defines  a  logical  aesh  nuaber,  KREGn  and/or 
LRBGn,  corresponding  respecitvely  to  a  physical 
distance,  XREGn  and/or  YREGn.  For  exaaple,  if 
XREG1  =  1.0  and  XMAX  =  2.0,  KRBG  =  21  and  KMAX  = 
61,  then  DXl  =  0.05  and  DX2  =  0.025. 


The  nuaber  of  curve  end  points  describing  the 
physical  boundary  of  the  region.  This  is  the 
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nuaber  of  data  iteas  used  to  describe  the 
boundary. 


* IPRINT=  * 
(0) 


If  IPRINT  =  1,  then  a  special  diagnostic  printout 
ia  provided  in  the  logical  path  finding  routine. 


Input  Variables  Describing  the  Physical  Coordinates  of  the 


Boundary 


Bach  data  line  gives  the  coordinate  of  an  end  point  of  a  curve  on 


the  physical  boundary  and  specifies  the  type  of  curve  to  be  drawn 


froa  the  previous  point  to  the  current  point.  Note  the  first 


data  line  is  the  starting  point  of  the  boundary. 


Coordinates  for  data  points  nay  be  given  as  (x,y),  (x’,y’),  (r, 


theta)  or  (r'f  theta')  and  aay  be  freely  interaixed  (see  Fig. 


2.1).  When  defining  an  arc  of  a  circle,  the  center  of  the  circle 


(X0,  Y0)  auat  be  given  and  the  end  point  of  the  arc  is  deterained 


relative  to  (X0,Y0).  An  exanple  of  an  arc  is  illustrated  in  Fig. 


2.2. 


The  following  is  a  reference  list  of  input  variables  that 


describe  the  physical  coordinates  of  the  boundary: 


NT,  XO,  YO,  X,  Y,  R,  THETA,  ICATH 


$ 
I® 


ft® 

'm 


k  complete  description  of  each  variable  is  given  below: 


*  MT=  * 

(1) 


*X0=* 
'  Y0=  * 
(0.0) 


*X=’ 
*  Y*  * 


’  R=  ’ 

*  THBTA=  * 


*  NEW=  * 

(0) 


’ ICATH=  * 

(0) 


Specifies  the  type  of  curve  to  be  drawn  froa  the 
previous  curve  endpoint  to  the  current  point. 

NT  =  1,  straight  line 

NT  =  2,  an  arc  of  a  circle, ( X0, Y0)  Bust  be  at  the 
center  of  the  circle. 

NT  =  3,  The  hyperbola  x’#  y’  =  (r  **  2)/2.0,  branch 
in  the  first  quadrant.  The  input  variable  HR" 
■ust  be  set  equal  to  the  r  in  this  equation. 


The  displacement  of  the  (x’,y’)  coordinate  axes, 
(0.0)  and  the  center  of  the  circle  if  "NT"  =  2. 


The  (X’,Y*)  coordinates  of  this  curve’s  end  point. 
If  NT=2,  these  values  are  on  the  new  (X’,Y’)  axis. 


The  (R*, theta’)  coordinates  of  this  curve  end 
point . 


NEW  =  1,  the  points  on  the  logical  path  for  this 
segment  are  not  allowed  to  coincide  with  the 
logical  boundary  points  of  any  previous  region. 


NEW  =  -1,  the  points  on  the  logical  path  for  this 
segment  are  not  allowed  to  coincide  with  the 
logical  boundary  points  of  any  previouys  region, 
except  for  the  starting  and  ending  points. 


Specifies  if  this  node  is  to  be  an  emitting  node 
or  not. 


=0  Non-emitting  node 
=1  Emitting  node 


2.6  Boundary  Condition  Data  Set-Up  by  AUTOMBSH 


The  boundary  condition  indicators,  C0N(21)  through  CON(24)  = 
NBSUP,  NBSLO,  NBSRT,  NBSLF,  are  set  by  AUTOMBSH  in  the 
LATTICB.DAT  "problem-constant"  data  as  follows: 


POISSON  Problen: 

NBSUP  =  0 
NBSLO  =  1 
NBSRT  =  0 
NBSLF  =  0 


SUPERFISH  Problen: 

NBSUP  =  1 
NBSLO  =  0 
NBSRT  =  1 
NBSLF  =  1 


The  region  special  boundary  indicators,  =  C(6)  in  the  LATTICE.DAT 
data  file,  are  set  by  AUTOMBSH  to  the  following  default  values: 


POISSON  Problem: 

C (6 )  =  1, 

condition) . 


for  all  regions  (Neumann 


SUPERFISH  Problem: 

C(6)  =  (1)  for  all  regions 

NOTE:  These  default  values  for  C(6)  may  be  over-written  by 

inputing  "IBOUND"  in  the  AUTOMESH  region  input  data.  See  Sec. 
3.6  of  LATTICE  for  a  description  of  the  special  region  boundary 
indicator. 
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2.7  Sample  Input  Data  for  AUTOMBSH.DAT  for  the  Problem  ’’MODGUN'* 

The  first  step  in  using  AUTOMESH  is  to  determine  the  coordinates 
of  various  corner  points  of  the  boundary  in  the  z-r  plane.  The 
dimensional  diagram  of  the  electrodes  for  the  Pierce  gun,  MODGUN, 
is  shown  in  Fig.  2.3.  The  boundary  of  MODGUN  with  it’s  corner 
points  is  shown  in  Fig.  2.4.  The  next  step  is  to  choose  a 
triangle  size  for  the  mesh.  This  is  accomplished  by  specifying 
the  input  variables  DX  and  DY.  Many  tines  it  is  only  necessary 
to  specify  DX  and  use  the  default  value  of  DY,  but  in  this 
example,  MODGUN  has  some  corners  which  require  finer  detail  than 
that  given  by  the  default  value  of  DY. 

Choosing  a  proper  DX  and/or  DY  is  important  so  that  AUTOMESH  can 
construct  a  consistent  logical  mesh.  If  DX  and/or  DY  is  too 
large,  AUTOMESH  will  not  be  able  to  find  a  logical  path  along  the 
boundary.  What  usually  happens  when  a  poor  DX  and  DY  has  been 
chosen  is  that  AUTOMESH  will  assign  the  same  logical  coordinate 
to  two  different  physical  coordinates  which  results  in  the 
program  being  aborted.  When  choosing  values  of  DX  and  DY,  a  good 
rule  of  thumb  is  not  to  use  a  value  much  larger  than  the 
smallest  difference  in  "X"  and  "Y"  distances  between  boundary 
corners . 


Once  DX  and  DY  has  been  chosen,  the  data  file  AUTOMESH.DAT  is 
created  by  the  user.  AUTOMESH.DAT  contains  the  boundary 
information  that  AUTOMESH  needs  to  construct  the  logical 
coordinates  of  the  boundary.  The  data  file  for  the  MODGUN 
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problem  is  shown  in  Fig.  2.5.  Line  nuaber  one  is  the  title  line, 
lines  two  through  four  are  the  input  variables  discussed  in  Sec. 
2.4  and  lines  five  through  forty-two  are  the  physical  coordinates 
discussed  in  Sec.  2.5. 

AUTOMBSH  creates  the  data  file  LATTICB.DAT  which  is  used  by  the 
progran  LATTICE  to  construct  the  aesh.  The  LATTICE. OAT  file  for 
the  problea  MODGUN  is  shown  in  Fig.  2.6  and  the  logical  diagraa 
constructed  by  AUTOMESH  is  shown  in  Fig.  2.7. 
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MOD -I  GUN 

I 

NRBG= *  1,  * DX= '  5.9B-2,  *DY=‘  4.33B-2,  ,XMIN=I  0.0, 
XMAX=*  2.8,  * YMIN  = '  0.0,  * YMAX  =  *  1.22,  *CONV=*  0.0254 
MPOINT=  *  18 

l 

X=*  0.0,  * Y= ’  0.0,  * ICATH= '  1 

X=*  2.7559,  * Y=  *  0.0 

X=*  2  7559,  * Y= ’  0.44600 

X=*  1.10800,  ' Y=  *  0.44600 

X=*  1.00800,  * Y* *  0.50150 

X=’  1.00800,  *  Y=  *  0.53850 

X=*  1.55500,  ' Y=  *  0.79700 

X=’  1.55500,  ’ Y=  *  1.20000 

X=’  1.03790,  ’ Y= ’  1.20000 

X=’  1.03790,  *  Y=  *  0.93500 

X=’  1.01400,  ’ Y=  *  0.89170 

X=’  0.96400,  *  Y=  *  0.89170 

X= ’  0.93890,  * Y=  *  0.93500 

X=*  0.93890,  * Y=  *  1.12900 

X=’  0.57290,  ’ Y= ’  1.12900 

X=’  0.28890,  *  Y= '  0.81300 

X=’  0.271058,  ’ Y= '  0.8020,  ’ICATH=’  1 

NT= ’  2,  ’X0=’  1.322,  ’Y0=*  0.0, ’X=’  -1.322,  *Y=*  0.0 
ICATH=  *  1 


.5  AUTOMESH . DAT  data  file  for  the  MODGUN  problen. 


HI*  LATTICE  (The  Hesh _ Generator) 

The  program  LATTICE  generatea  the  non-uniform  triangulated  meah 
from  input  data  daacribing  the  boundary.  Input  data  ia  read  by 
the  format-free  input  routine  FREE  which  ia  described  in  Sec. 
1.4.  The  terminology  used  throughout  these  notes  is  as  follows: 
an  "air"  region  means  a  vacuum  region  in  which  the  relative 
permeability  equals  1.0,  and  an  "iron"  region  is  a  region  which 
contains  material  for  which  the  relative  permeability  is  a 
constant  but  not  equal  to  1.0,  or  defined  by  a  linear  function  or 
a  non-linear  table. 

U _ Icpat  Data  for  LATTICE 

There  are  five  groups  of  input  data  for  LATTICE: 

1)  The  first  data  line  (card)  can  have  anything  in  columns 
2  through  80  and  these  characters  are  stored  for  run 
identification  in  the  printed  output.  If  column  1  is 
blank,  then  this  data  is  for  a  POISSON  problem.  If  the 
column  is  non-blank,  then  this  data  is  for  a  SUPERFISH 
problem.  This  is  shown  in  Fig.  3.1a. 

2)  The  second  data  line  (or  lines)  give  values  for  the 
"problem-constants''  (CON’s).  These  are  general 
parameters  and  options  specifying,  for  example,  the 
number  of  regions,  the  boundary  conditions,  the  type  of 
magnetic  symmetry,  etc.  An  input  value  for  C0N(2)  is 
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EXAMPLE  LATT1CE.DAT  RLE 


U 


•2  2  *21 

1  1  0 

0  *9  1.0000 

•138  4  SKIp 

EE 

0.0000 

0.0000 

0 

0  REGION 

] 

rrr 

0.0000 

0.0000 

B 

1 

10  1 

1.0000 

0.0000 

B 

10  10 

1.0000 

1.0000 

B 

1  10 

0.0000 

1.0000 

B 

1  1 

L 

0.0000 

0.0000 

B 

COUNT 

- 1 


B(0  ZB(I) 


Fig.  3.1 


The  five  data  groups  of  the  LATTICE  input  data 
file. 

aj  Title  line. 

bj  Values  of  the  "problem-constants",  CON’s. 

c)  "region-constants",  C(l)  array. 

d)  Logical  and  physical  boundary  coordinates, 
e;  Axial  magnetic  field  data  (optional). 
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is  only  required,  the  others  all  having  default  values. 
The  foriat  of  this  line  is  illustrated  in  Fig.  3.1b 

and  the  input  list  is  terninated  with  the  character 
"S". 

Data  for  up  to  six  properties  of  the  region  nay  be 
specified  in  the  third  group  of  input  data.  In  the 
progran,  values  are  stored  in  the  "C"  array  and  these 
"region-constants"  ("C")  specify  the  following 
infornation: 

C(l)  The  region  nunber  (integer).  The  naxinun  nunber 
of  regions  is  30  and  they  need  not  be  consecutive. 

C(2)  The  naterial  code  for  the  region  (integer). 

=0  All  points  inside  this  region  are  onitted 
fron  the  problen. 

=1  air  (nu=l) 

=2  iron  with  internal  peraeability  table 

=3  iron  with  input  peraeability  table  no.  1 

=4  iron  with  input  permeability  table  no.  2 

=5  iron  with  input  peraeability  table  no.  3 

C(3)  The  fixed  potential  value  (volts)  on  the  boundary 

of  the  region  for  electrostatic  probleas,  and  for 

aagnetostat ic  probleas  it  is  the  total  current 

(saps)  in  the  region. 

=0.0  default  value  (real  nuaber) 
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C(4)  The  charge  density  (Couloabs/length**2 )  of  the 
region  for  electrostatic  problems  and  the  current 
density  (aaps/length**2)  of  the  region  for 


■agnetostatic  problems. 


=0.0  default  value  (real  nuaber) 


C(5)  Triangle  type  to  be  generated  in  the  region 
( integer) . 


=0  equal  weight  (default) 
= 1  equilateral 
=2  right 


C(6)  The  region  boundary  indicator  (integer).  See 
section  3.6  for  a  detailed  description  of  the 
boundary  conditions. 


=-l  Fixed  potential  Oirichlet  boundary. 

=  0  Magnetic  Oirichlet  condition  on  the  boundary 
of  this  region. 

*  1  Neuaann  condition  on  the  boundary  of  this 
region . 


The  third  group  of  data,  "region-constants",  is  shown  in  Fig. 
3.1c. 

4)  The  fourth  group  of  input  data  consists  of  the  following 
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three  iteas: 

1.  logical  coordinates  (K,L)  of  the  boundary. 

2.  physical  location  (x,y)  of  the  logical  coordinates  of 
the  boundary. 

3.  node  type 

The  format  of  this  input  data  in  LATTICE.DAT  appears  as: 

K( 1 )  L(l)  X(l)  Y ( 1 )  B 
K(2)  L ( 2 )  X ( 2 )  Y ( 2 )  K 
•  •  •  •  • 

•  •  •  •  • 

K(n)  L(n)  X(n)  Y(n)  K  Count 

The  input  list  is  terminated  with  the  character  "C "  to  count 
the  number  of  input  values.  Comments  may  follow  this  "C" 
since  everything  else  on  this  line  will  be  ignored.  The 
fourth  data  group  is  shown  in  Fig.  3. Id. 

Bach  logical  coordinate  must  include  it’s  node  type. 
Presently  there  are  three  node  types  which  are:  boundary, 
guard  and  cathode  points. 

The  fifth  group  consists  of  axial  magnetic  field  data 

(optional).  LATTICE  accepts  values  of  B  measured 

2 

experimentally  along  the  z-axis  from  which  off  axis  values 

of  B  and  B  are  calculated, 
z  r 

If  the  variable  NBPTS  (C0N(138))  is  greater  than  zero,  then 
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the  following  FORTRAN  format  free  READ  statement  is 
executed: 


READ (91,*) (B(I) ,ZB(I) ,I=1,NBPTS) 


where : 


B ( I )  -  value  of  the  magnetic  field  on  the  z-axis  (in 

Tesla) 


ZB(I)  -  z-axis  location  of  the  magnetic  field  B(I). 


The  fifth  group  of  data,  the  axial  magnetic  fields,  is  shown  in 
Fig.  3.1e.  B ( I )  and  ZB(I)  is  located  on  the  data  line  immediately 
following  the  last  region  and  there  should  be  exactly  NBPTS 
pieces  of  data  for  both  B(I)  and  ZB(I).  Also  note  that  the 
magnetic  field  data  should  extend  beyond  the  boundary  by  a  small 
amount . 
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Data  groups  (3)  and  (4)  are  repeated  for  each  region.  The  first 
reg  ion  defines  the  entire  problem  mesh  and  must  contain  the 


maximum  (K,L)  values  of  the  mesh.  Data  for  the  succeeding 
regions  redefines,  or  overwrites  the  region  constants  for  all 
triangles  belonging  to  this  new  region.  Data  for  each  new  region 
overwrites  previously  defined  values  in  the  same  way.  This 
feature  is  known  as  "successive  region  data  overwriting".  For 
example,  suppose  that  one  wishes  to  specify  the  cathode  at  one 
voltage  and  the  focus  electrode  at  another.  Two  additional 
regions  must  now  be  added  that  specify  the  cathode  at  and  the 
focus  electrode  at  . 

Usually  each  region  is  closed,  i.e.,  the  data  for  the  first  and 
last  boundary  points  are  identical.  It  is  also  possible  to 
specify  data  for  a  "point”  region  or  a  "line"  region.  The 
purpose  of  this  would  be  to  define  the  physical  coordinates  of 
specific  points  such  as  fixed  potential  values  along  a  portion  of 
a  boundary.  For  point  and  line  regions,  the  input  values  for 
C(2)  and  C(5)  in  the  "region-constant"  data  are  arbitrary,  since 
this  information  is  used  only  for  closed  regions. 

3.2  Specifying  Logical  Coordinates  In  Logical  Space 

For  two  points  in  logic  space,  it  is  necessary  that  they  be  on 
the  same  "logical  lines”  in  logic  space.  These  four  lines  shown 
in  Fig.  3.2  are: 
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Fig.  3.2 


Primary  "Logic  Lines"  on  a  Equilateral 
Triangle  Logic  Mesh 


w 


1. 


vertical 


2.  horizontal 

3.  slant  right 

4.  slant  left 

It  is  only  necessary  to  specify  the  end  points  of  the  logical 
line,  interaediate  points  will  then  be  linearly  interpolated  for 
(x.y) . 


3.3  Mesh  Generation  Resulting  In  Negative  Areas 


■■mwSsflffli nmmmn 
WWrATATATATATATATATA'V 
I  wwATiViTiViVin  I 
lBgwro»yiMrAwATial 


(a)  A  mesh  with  negative  area  triangles 


(b)  Corrected  mesh 


Fig.  3.3 


(a)  An  example  of  a  square  mesh 
that  has  some  negative  area 
triangles,  (b)  Corrected  square 
mesh  by  specifying  the  lower  left 
hand  corner  at  (0.0, 0.0) 
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1.  physical  coordinate  (x,y)  data  error 


2.  too  many  aesb  points  defined  in  a  certain  area  of  the 


■eah. 


For  the  exanple  of  Fig.  3.3a,  the  cause  of  the  error  was  case  *1. 


Instead  of  the  lower  left  corner  being  specified  at  (0.0,  0.0), 


it  was  specified  in  the  center  of  the  nesh,  which  resulted  in  the 


lesh  folding  over  upon  itself.  An  exanple  of  case  #2  will  be 


illustrated  in  the  exanple  of  MODGUN  in  Sec.  3.8. 


The  following  procedure  can  be  used  to  correct  the  problen  of 


negative  areas: 


Deternine  the  approxinate  physical  location  (x,y)  of 


the  negative  area  triangle.  This  can  be  deternined 


fron  the  knowledge  of  the  (K,L)  coordinate  and  the 


LATTICE . DAT  data  file. 


2.  Using  the  plotting  progran,  zoos  in  on  the  troubled 


area  of  the  nesh. 


3.  After  examining  the  problen,  deternine  what  changes 


nust  be  nade  to  LATTICE.DAT  to  correct  the  problen. 


As  in  the  exanple  of  Fig.  3.3a,  noving  the  lower  left  hand  corner 


fron  the  center  of  the  nesh  to  (0.0,  0.0)  corrected  the  problen. 


This  is  shown  in  Fig.  3.3b. 
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The  phyeicel  location  (x,y)  of  each  aeah  point  (or  node)  within 


the  boundary  of  a  problem  is  determined  numerically  by  a 


relaxational  process.  In  some  instances,  the  user  may  want  to  fix 


certain  mesh  points  at  a  predetermined  location  (x,y)  in  the 


mesh.  The  "guard"  node  is  used  to  fix  the  (x,y)  coordinate  of  a 


node. 


The  "guard"  mesh  is  a  separate  region  in  the  LATTICE.DAT  file. 


The  "guard”  region  should  be  defined  after  the  entire  problea 


meah  of  the  firat  region  has  been  defined.  For  "guard"  regions, 


the  input  values  for  C(2)-C(5)  in  the  "region-constant"  data  are 


arbitrary  and  C(6)  can  either  be  set  to  C(6)=0  or  1,  but  not  -1. 


An  example  illustrating  the  uae  of  "guard”  nodes  is  shown  in  Fig. 


3.4.  The  "guard"  mesh  is  set  up  on  the  X  =  0.1  line. 


mentioned  in  Sec.  3.2,  it  is  only  necessary  to  specify  the  end 


points  of  the  logical  line,  intermediate  points  are  linearly 


interpolated  for  (x,y). 


-  40  - 


a)  Square  mesh  without  "guard"  line. 
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b)  Square  mesh  with  "guard"  line. 


Fig.  3.4  Example  of  a  square  mesh  with  a  "guard"  line 
at  x  =  0.1.  a)  IATTICE.DAT  data  file  and 
relaxed  mesh  without  the  "guard"  line, 
b)  LATTICE.DAT  data  file  and  relaxed  mesh  with 
the  "guard"  line. 


The  emitting  surface  conaiata  of  a  aet  of  eleaentary  planar 


diodes.  Near  this  surface  an  approximately  rectangular  meah  ia 
generated  as  shown  in  Fig.  3.5.  The  spacing  "z"  of  the 
rectangular  nesti  is  determined  by 

z  =  (a+b)  /  ( FRTtNAC ) 

where  z,  a  and  b  are  defined  in  Fig.  3.5.  The  program  parameter 
F8T  (CON(137)>  controls  the  nearness  of  the  nodes,  and  the 
program  parameter  NAC  (CON(126))  is  the  number  of  nodes  in  front 
of  the  emitting  surface. 

It  has  been  determined  empirically  that  for  curved  emitting 
surfaces  being  approximated  by  a  aet  of  planar  diodes,  two 
criteria  should  simultaneously  be  : 

1.  "The  change  in  the  angle  of  the  line  of  nodes  in  the 
rectangular  mesh  from  cathode  point  to  cathode  point 
ahould  be  less  than  4  degrees". 

2.  The  program  parameter  FRT  which  controla  the  nearness 
of  the  nodes  should  be  greater  than  8. 

The  user  must  be  careful  that  these  two  conditions  are  satisfied 
since  the  potentials  deviate  substantially  from  planar  diode 
theory  when  they  are  not  close  to  the  emitting  surface. 
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The  nodes  "N"  of  Pig.  3.5  are  the  slitting  points  of  the  cathode. 
They  are  designated  as  "I"  nodes  in  the  LATTICB.DAT  data  file. 

3.5.1  Horizontal  and  Vertical  "Lozic-Line"  Cathodes 

The  prograa  can  generate  two  basic  types  of  cathode  sesh  in 
"logic"  space: 


1.  vertical 

2.  horizontal 


The  "horizontal"  cathode  aust  be  specified  on  the  sase  horizontal 
logic  line,  i.e.,  all  cathode  nodes  sust  be  on  the  ease  "L"  line 
in  logic  space.  The  NAC  nodes  associated  with  the  cathode  nodes 
are  assigned  froa  bottoa  to  top  on  the  sane  T  line  as  shown  in 
Fig.  3.6a.  In  cylindrical  coordinates,  a  horizontal  cathode  on 
the  "z"  axis  is  not  allowed  since  this  situation  is  not 
physically  possible. 

The  rules  for  the  "vertical"  cathode  are  less  stringent  than 
those  of  the  "horizontal”  cathode.  The  "vertical"  cathode  »ay  be 
constructed  froa  a  vertical,  slant  right,  slant  left  or  a 
coabination  of  these  priaary  "logic  lines”  (see  Fig.  3.2).  The 
NAC  nodes  associated  with  the  cathode  nodes  are  assigned  from 
left  to  right  on  the  ease  ”L"  line  as  shown  in  Figs.  3.6b,  c  and 
d. 

When  constructing  a  "vertical  logic-line"  cathode,  it  is 
iaportant  that  each  eaitting  node,  "K",  has  NAC  nodes  free  to  the 

-  44  - 


!»s 


m 

m 


Emitting  NodM  ■  i 
Cothode  M«ah  Nodes  -  O 


w^^AmwAWArAW^AFAW^AWAri 

Ltatatatatatatatatatatatatatatatatatatatai 

r^ATATATATArATATATArATArATATATATATATATAfi 
OTAWAWAYaVATaWATAWATATATaWAW 
AT^^^ArAWAfATATAW^ATATATATArArAfi 
TATA  iwmw  atatatatatatatatatay 
riMNIA%nnnn^ArATArATATATATATATi^H 
TAW  rlV#V.VoW  AWATaYaYAYAYATOS 
ATAMktAtAtAtAtArA^TATATATATATATATATAU 
ATATA  ATATaTATaWA^A 
rATATATATATAYArATArATATATATATATATATATATAri 


Logical  Path  of 
tho  Boundary 


Cathode  M«h 
Area 


Fig.  3.6  a)  Logical  space  diagram  and  relaxed  triangular  mesh 
of  the  horizontal  cathode. 
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Fig.  3.6  b)  Logical  space  diagram  and  relaxed  triangular  mesh 
of  the  vertical  cathode. 
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Fig.  3.6  c)  Logical  space  diagram  and  relaxed  triangular  mesh  of  the 
vertical  slant  left  cathode. 
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Fig.  3.6  d)  Logical  space  diagram  and  relaxed  triangular  mesh  of  the 
vertical  slant  right  cahtode. 
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right  of  it,  i.e.,  no  two  emitting  nodes  may  be  on  the  sane  "L” 
logic  line  (or,  no  two  emitting  nodes  may  be  on  the  same  "K"  line 
for  a  "horizontal  logic-line"  cathode).  The  program  LATTICE 
does  not  check  for  this  error,  and  as  a  result,  negative  area 
mesh  will  be  generated.  An  example  of  a  "slanted  logic-line" 
cathode  that  shares  the  same  logic  "L"  value  for  two  emitting 
points  is  shown  in  Fig.  3.7a,  and  the  corrected  mesh  is  shown  in 
Fig.  3.7b. 
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The  user  aust  be  cautioned  +  fcat  AUTOMESH  will  not  construct  a 
cathode  " logic-1 ine”  according  to  the  rules  given  in  Sec.  3.5.1. 
As  discussed  in  Chapter  2,  AUTOMESH  finds  the  path  in  the  logical 
■eah  which  best  fits  the  physical  boundary.  This  logical  path 
found  by  AUTOMESH  will  not  necessarily  cohere  to  the  rules  of 
Sec .  3.5.1. 

Probably  the  best  way  to  construct  a  consistent  cathode  mesh  from 
the  output  of  AUTOMESH  is  to  put  all  the  cathode  ”K"  nodes  on  one 
primary  logic  line.  As  for  the  example  MODGUN  of  Sec.  2.7,  the 
cathode  could  be  changed  from  the  "zig-zag”  vertical  logical 
line,  to  the  K=1  line.  This  is  shown  in  Fig.  3.11  of  Sec.  3.8. 

Another  area  of  concern  on  the  cathode  is  the  last  emitting 
point.  Many  times  AUTOMESH  will  put  a  boundary  point  on  the  same 
logical  line  "l"  as  the  last  emitting  point.  As  a  result,  some 
boundary  points  are  lost  when  the  cathode  mesh  is  generated  by 
LATTICE  because  the  last  emitting  point  did  not  have  NAC  free 
nodes  to  the  right  of  it.  An  example  of  this  type  of  problem  is 
shown  in  Fig.  3.8a.  This  problem  was  corrected  by  changing  the 
boundary  point  from  logical  coordinate  (20,2)  to  (21,1).  The 
results  of  the  correction  is  shown  in  Fig.  3.8b  (Note:  Fig.  3.8b 
includes  a  "guard”  mesh). 
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b)  Corrected  meah  by  moving  the  boundary  point  above  the  lost  emitting 
node. 


Fig.  3.8  Modifying  incorrect  boundary  point  generation  of  AUTOMESH 
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There  are  two  ways  that  boundary  conditions  are  specified  in 


LATTICE.  These  two  Methods,  specifying  CON ( 2 1 ) -CON ( 24 )  and  C(6) 
of  the  "region-constants” ,  are  described  below. 


Sec.  3.6.1  Specifying  CON ( 2 1 ) -CON( 24 ) 

The  four  CON’s,  C0N(21),  CON(22),  CON(23),  and  CON(24), 

determine  the  boundary  conditions  on  the  extreme  rectangular 
logical  boundaries  of  the  problem  regardless  of  the  value  of  C(6) 
for  a  region  part  of  whose  boundary  coincides  with  an  extreme 
rectangular  logical  boundary.  The  extreme  logical  boundaries  are: 
K= 1 ,  and  K  =  KMAX ,  and  L  =  1  and  L  =  LMAX .  (Note:  if  C(6)  =  -1, 
then  this  is  not  the  case).  See  Sec.  3.7.1b  for  the  description  of 
these  CON’s. 

Sec.  3.6.2  "Region-Constant"  C(6) 

C(6),  of  the  "region-constant"  C  array  described  in  Sec.  3.1,  is 
used  for  two  purposes: 


1)  to  indicate  the  boundary  condition  on  a  boundary  of  the 

problem  which  does  not  coincide  with  the  extreme 
logical  boundaries  (K=l,  1=1,  K  =  KMAX,  L=  LMAX )  of  the 

problem. 

2)  to  indicate  special  fixed  potential  points. 
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1)  When  used  for  purpose  1,  C(6)  acts  as  either  a  Neumann  or 

Dirichlet  boundary  indicator: 


C(6)=0  Dirichlet  condition  on  the  boundary  of  the 

particular  region  (  sagnetic  Dirichlet 
condition  ). 

C(6)=l  Neusann  condition  on  the  boundary  of  the 

particular  region. 

The  default  values  for  C(6),  which  are  used  unless  changed  by 
input  data  are: 

C(6)=l  for  all  regions. 

2)  When  C(6)  is  used  for  purpose  (2),  it  is  set  to  -1  and  the 
fixed  potential  value  is  C(3)  of  the  "region-constant"  data. 

When  using  feature  (2),  the  region  should  be  either  a  line  or 
point  region  that  overwrites  a  portion  of  the  boundary.  C(6) 
should  never  be  set  to  -1  for  the  first  region,  which  describes 
the  whole  boundary. 

When  defining  the  boundary  of  a  problee  (  first  region  ),  C(6) 

should  be  set  to  one,  and  all  Dirichlet  boundaries  should  be 
added  in  succeding  regions. 

Sec.  3.8  gives  an  exaaple  of  using  both  boundary  indicators  and 
the  use  of  line  regions  to  define  fixed  potentials  along  the 
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Sec.  3.7  "Prob lea-Constant"  Arra 

The  ’’CON"  array,  ia  an  array  containing  integera  and  real  nuabers 
which  are  general  paraaetera  and  options  specifying  the  control 
of  the  prograa.  These  values  are  specified  or  changed  from  their 
default  values  in  the  second  group  of  input  data  for  LATTICB  (see 
Sec .  3.1). 


t 


k 


Sec.  3.7.1  describes  the  CON(I)'s  which  controls  the  aesh 
generation  of  LATTICB.  Sec.  3.7.2  describes  the  CON(I)’s  that 
are  associated  with  the  electron-flow  (gun)  problem. 

Sec.  3.7.1  Mesh  Generation 

This  section  describes  the  ’’problem-constants'*  that  can  only  be 
specified  in  the  LATTICB  program.  These  CON’s  are  general 
parameters  and  options  specifying  the  control  of  the  aesh 
generation  in  LATTICB.  Default  values  are  given  in  parentheses. 


3.7.1a  Parameters  defining  the  basic  problem: 


CON ( 2 ) =NRBG 


Number  of  regions  in  the 
LATTICE  data  file. 


CON ( 9 ) =CONV 
(1.0) 


Conversion  factor  for  the 
units  in  the  LATTICB  data 
file. 

CONV  =  1.0  for  meters 
CONV  =  0.0254  for  inches 
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3.7.1b  Boundary  Condition* 


C0N(21 ) =NBSUP 

(0) 

CON(22)=NBSLO 

(0) 

CON (23)=NBSRT 
(0) 

CON (24)=NBSLF 

(0) 


Indicators  for 
boundary  conditi 
UPper,  LOwer,  Ri 
logical  bounder 
problem  region,  e 
Neuaann  boundary 
logical  problea 
NBSUP=0 ,  Dirichle 
the  upper  logi 
boundary. 


the  Neuaann 
ons  on  the 
ghT,  and  LeFt 
ies  of  the 
.g.,  NBSUP=1 , 
on  the  upper 
boundary , 
t  boundary  on 
cal  problea 


3.7.1c  Printed  Output  Options: 


CON(32)=IPRINT 

(0) 


IPRINT=- 1 ,  print  (x,y) 
coordinates  of  aesh  points. 
IPRINT=0,  nothing 


3.7. Id  Current  Paraaeter  and  Options: 


CON ( 70 ) = ICAL  ICAL= 1 ,  use  angle  foraula  for 

calculating  the  current 

associated  with  a  point  (used 
when  accurate  fields  near  coil 
boundaries  are  needed). 

ICAL=0,  use  noraal  area 
foraula. 


3.7. le  Over  Relaxation  Factors 


C0N(79)=RH0XY 

(1.6) 


The  starting  over-relaxation 
factor  for  the  generation  of 
the  aesh. 


CON ( 81 ) =N0TE 

(1) 


NOTB=0,  point  relaxation  order 
is  air  points,  interface 
points,  then  iron  points  (this 
point  order  aust  be  used  for 
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PANDIRA) 


N0TB=1,  point  relaxation  order 
ia  air  and  interface  points, 
then  iron  points. 


3.7. If  Convergence  Criteria 


CON ( 84 ) =BPSO 
(1.0B-5) 


The  convergence  criterion  for 
the  aesh  generation. 


3.7. lg  Conforaal  Transf oraat ion  Factors 


CON ( 37 ) =MAP 

(1) 


CON( 123) =TNBGC 
CON( 124 ) =TPOSC 


W=Z**MAP/ (MAPARZERO** (MAP- 1 ) ) 
is  the  conforaal 

transf oraat ion  used,  if  any. 

Prograa  LATTICE  autoaat ical ly 
transforaa  the  current 

density.  Prograa  POISSON 

transforas  the  peraeability 
and  calculates  fields  in  both 
the  original  and  the 

transforaed  geoaetrics. 


The  total  negative  and 
positive  current  in  the 
transforaed  geoaetry  (this 
equals  the  total  negative  and 
positive  current  in  the 
original  geoaetry). 


CON( 125) =RZBRO 
(1.0) 


The  scaling  factor  if  the 
transforaation 

*f=Z*»MAP/MAP*RZBRO**(MAP-l)  ) 
is  used.  Nornally 

RZBR0=ABS (ZO)  which  equals  the 
aperture  radius. 
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Sec.  3.7.2  Space-Charge  Flow  Parameter* 


This  section  describes  CON’s  (125)-(150)  which  are  associated 


with  the  electron-gun  problem. 


C0N(126)=NAC 

(3) 


CON( 127) =SSSF 

(0.2) 


C0N( 128) =CGAMMA 
(0.2) 


CON (129)=FJLIM 
(  IB  10  ) 


C0N( 130)=CHILDK 
(2.334E-6) 


C0N( 131 ) =V0LTCA 

(0.0) 


C0N( 132 ) = IGUN 


Nuaber  of  nodes  in  the 
rectangular  aeah  in  front  of 
the  cathode. 


Space-charge  suppression  factor 
(used  when  CBAR<=0) 


CBAR  modification  parameter, 
or  cathode  loading  convergence 
parameter . 


Maximum  current  density  from 
temperature  limited  operation 
considerations  (asps/i  ) 


Child’s 

constant 


planar  diode 


Cathode  to  Anode  voltage  (this 
is  automatically  determined  in 
POISSON) . 


Control  flag  indicating 

electron-gun  problem. 


IGUN=0, 

problem. 


electron-gun 


IGUN=1,  electron  gun  problem 
(LATTICE  will  automatically 
determine  this  by  the  presence 
of  cathode  "K"  nodes). 


m 


w- 


8 


bisssi 


(MJurwiuHiwtviHnjumjiiuvuouMJuriKnuMnjuuuvui 


CON(133)~STA 

(1.759811) 


CON( 135) =MAXDT 
(400) 


CON( 136) = NODES 


CON(137)=FRT 

(8.0) 


CON( 138) =NBPTS 

(0) 


CON( 139) =FRDT 

(0.12) 


CON ( 140 ) =MAXCYC 
(20) 


CON (141) = EC BAR 
(5. OB-3) 


CON( 142) =PQUBSS 

(0) 


CON ( 143 ) =MPRI NT 

(1) 


CON ( 144 ) =DCBAR 


Blectron  charge  to  Bass  ratio. 


Maxiaua  nuaber  of  tiae  steps 
allowed  in  the  tracing  of  each 
ray. 


Nuaber  of  eaitting  nodes  of 
the  cathode  (deterained  in 
POISSON) . 


Paraaeter  which  controls  the 
nearness  of  the  rectangular 
aesh  in  front  of  the  cathode. 


Nuaber  of  aagnetic 
points  to  be  inputed 
NBPTS=0,  no  aagnetic 
data  is  read  in). 


field 

(for 

field 


A  aultiplication  constant  in 
which  the  time  step  say  be 
fractionally  adjusted. 


Maxiaua  nuaber  of  aajor 
current  cycles  to  be  allowed 
to  run. 


CBAR  convergence  tolerance 


( anps/vo 


perveance 


Printing  options  for  the 
results  of  POISSON. (see  4.4.2) 


The  change  in 
aajor  current 


CBAR  between 
cycles .  This 
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quantity  is  calculated  in 
POISSON,  and  is  compared  with 
ECBAR  to  determine 
convergence. 


CON ( 145 ) =CURREN 

(0) 


The 
from 
cathode . 


amount  of  current  emitted 
a  temperature  limited 


CON (146)=INCPLT 
(1) 


Time  step  increment  for 
plotting  and  printing 
trajectories . 


CON ( 147 ) = ISP 
(8) 


Spiral  scan  limit  of  the  block 
type  search  in  the  nearest 
node  routine. 
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Sec.  3.8  Sample  LATTICE  Data  File  For  The  Problem  MODGUN 


In  Sec.  2.7,  AUTOMESH  created  the  data  file  LATTICE.DAT  for  the 
problem  MODGUN  as  shown  in  Fig.  2.6.  A  minimum  of  three 
Modifications  have  to  be  nade  to  LATTICE.DAT  before  the  aesh  can 
be  generated  by  LATTICE: 

1)  Set  C0N(2)  and  CON( 21 ) -CON ( 24 )  to  their  correct  values. 

C0N(2),  which  is  the  number  of  regions,  is  set  to 
three.  The  three  regions  are:  problem  region,  cathode 
and  focus  electrode  region,  and  the  anode  region. 
C0N(21),  C0N(22)  and  CON(23)  are  set  to  one  since  the 

axis,  the  extreme  top  logical  L  line  (LMAX),  and  the 
extreme  right  logical  K  line  ( KMAX )  are  Neumann 
boundaries . 

2)  Change  the  "zig-zag"  vertical  logic  line  to  a  constant 
K=1  logic  line.  The  "zig-zag"  logic  line  would  have 
worked,  but  for  the  sake  of  an  example  in  constructing 
a  consistent  cathode  mesh,  this  modification  is 
illustrated. 

3)  Add  two  additional  regions: 

1)  cathode  and  focus  electrode  at  0  volts 

2)  anode  at  8000  volts 

Note  that  IB0UND=C(6)  is  set  to  -1  for  these  two 


The  aesh  was  generated  with  these  three  aodif ications,  and  a 
negative  area  aesh  occurred  at  the  top  edge  of  the  cathode  as 
shown  in  Fig.  3.9a.  The  folding  of  the  aesh  occurred  due  to  the 
close  proxiaity  of  the  corners  and  the  large  nodal  density  of  the 
cathode  aesh.  This  problea  was  corrected  by  using  a  "guard"  line 
as  shown  in  Fig.  3.9b. 

The  input  data  file  for  LATTICE,  including  the  "guard"  regions, 
is  shown  in  Fig.  3.10  and  the  logical  diagraa  is  shown  in  Fig. 

3.11.  The  aesh  generated  for  the  problem  MODGUN  is  shown  in  Fig. 

3.12. 

Note  the  addition  of  the  two  "guard"  regions  in  Fig.  3.11  and 
Fig  3.12.  These  "guard"  regions  were  added  to  illustrate  the 
use  of  a  "guard"  aesh  in  a  problea. 
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Fig.  3.10  Data  file  LATTICE.DAT  for  the 
problem  MODGUN. 
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IV.  POISSON  (Solution  of  the  Space-Charge  Flow  Problem) 

The  prograa  POISSON  solve*,  by  successive  point  over-relaxation, 
Poisson's  equation  in  two  dimensions.  This  portion  of  the  aanual 
describes  the  procedures  for  solving  the  space-charge  flow 
problea.  For  aagnetostat ic  problems,  the  user  should  refer  to 
the  original  aanual  for  instructions  to  solve  this  type  of 
problea. 

4.1  Input  Data  for  POISSON 

There  are  two  groups  of  input  data  for  the  space-charge  flow 
problea  and  three  groups  of  input  data  for  the  electrostatic 
problea: 

1)  The  first  group  of  input  data  is  an  integer  number, 


"NUH" ,  which 

is  the  nuaber 

of 

the  "dump”  to 

be  read  and 

processed  by 

POISSON. 

For 

space-charge 

flow  and 

electrostatic 

probleas. 

the 

duap  number 

0  (zero) 

written  by  LATTICB  is  to  be  used  (  see  Fig.  1.3  of  Sec. 
1.5  ).  Dump  numbers  greater  than  zero  are  written  by 
POISSON  and  duaps  greater  than  one  are  associated  with 
magnetostatic  probleas. 

2)  The  second  group  of  input  data  are  the  input  values  for 
the  "problem-constants” ,  CON’s.  These  input  values  can 
change  previously  defined  CON  values  read  froa  duap  0, 
or  specify  values  which  are  used  only  in  POISSON  and 
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were  not  input  in  the  LATTICE  problem-constant"  data. 
The  foraat  of  the  second  group  of  input  data  is  read  in 
by  the  foraat  free  input  routine  and  is  described  in 
Sec.  1.4.  This  foraat  is  exactly  the  saae  as  data 
group  (2)  of  LATTICB. 


The  third  group  of  data  is  a  flag  to  terainate  the 
POISSON  run  and  is  associated  with  the  electrostatic 
and  aagnetostat ic  problea.  The  last  data  line,  -1,  is 
typed  to  the  terainal  to  terainate  the  run. 


POISSON  "Problem-Constants" (CON's) 


This  section  describes  the  CON’s  which  controls  the  solution  of 
Poisson’s  equation.  Some  CON’s  apply  to  aagnetostat ic  problems, 
and  the  user  should  refer  to  the  original  annual . 


4.2.1  Description  of  the  CON’s 


The  following  sections  describe  the  CON’s  used  in  POISSON.  The 
default  values  of  the  CON’s  are  shown  in  parentheses. 


4.2.1a  Parameters  Defining  the  Basic  Problea 


CON ( 6 ) =MODE 
(-2) 


Code  for  the  peraeability  in  iron. 
Mode=-2,  au-infinite  in  iron 
Mode=-l,  au-finite  but  constant 
using 

gamma =1 . 0/mu= CON ( 10 ) = FI XGAM . 
Mode=0,  mu- f i n i t e- var i ab le ,  using 
the  internal  (B, gamma)  table 
or  a  table  input  by  the  user. 


C0N(7)*STACE 

(1.0) 


Stacking  (fill)  factor  in  iron 
ragiona . 


CON( 10)=FIXGAM 
(0.004) 


The  value  of  gaaaa  uaed  in  a  au- 
finite  but  conatant  aolution 
(CON (6) =M0DB=-1 ) .  Alao  the  value 
to  which  gaaaa  ia  initialized  for  a 
au— finite- vari ab 1 e  aolution. 


C0N( 18) SNPBRM 

(0) 


The  nunber  of  peraeability  tablea 
to  be  read  in  aa  data  by  prograa 
POISSON. 


CON( 19)=ICYLIN 

(1) 


Indicator  for  cylindrical  syaaetry 
problea. 

ICYLIN=1,  cylindrical  ayaaetry 
(r,*) 

ICYLIN=0,  two  diaenaional  (x,y) 


_ Boundary  Conditions 


CON ( 20 ) = INPUTA 

(0) 


The  nuaber  of  fixed  potential 
valuea  to  be  read  in  aa  data  by 
prograa  POISSON. 


CON(50) =IHDL 
(100000) 


During  the  POISSON  iteration,  the 
nuaber  of  cyclee  between  asking  a 
quasi-integral  Hdl  calculation 
around  the  Dirichlet  boundary. 
This  soaetiaes  speeds  the 
convergence,  particularly  for  non- 
syaaetrical  "H"  aagnets. 


4.2.1c  Liaita.  Printed  Output  and  "Duap"  Options 


CON ( 30 ) =MAXC YC 
( 100000) 
CON( 31 ) = IPRFQ 


Maxiaua  nuaber  of  iteration  cycles. 
The  cycle  print  frequency  during 


(0)  (0)  the  POISSON  iteration.  The 

default  value  of  0  indicates  that 
the  iteration  information  will  be 
printed  only  on  the  first  and  last 
cycles.  Input  values  of  IPRFQ  Bust 
be  an  integer  multiple  of 
I VERG=C0N ( 87 ) . 


C0N( 32 ) IPRINT  IPRINT=-1,  Print  (x,y)  coordinates 

(0)  of  the  mesh  points  (LATTICE  only). 

IPRINT=Q ,  nothing 

IPRINT=1,  Print  the  vector 
potential  array  when  CON ( 6 ) =M0DE =0 . 
IPRINT=2,  Print  the  ;Bi  in  iron 
triangles . 

IPRINT=4,  Print  the  B  ,B  )  in  iron 
triangles.  x  y 

Any  combination  of  these  options 
may  be  specified  by  the  sum  of 
"IPRINT”  values. 


CON( 35 ) =NODMP 
(0) 


Indicator  not  to  write  a  "dump"  at 
the  completion  of  a  POISSON 
solut ion . 

NODMP=0 ,  Write  "dump" 

N0DMP=1,  Do  not  write  "dump" 


4. 2. Id  Bdit  Options  (Fields  and  Gradients) 


CON ( 38 )  =  XORG 
(0.0) 

CON ( 39 ) = YORG 
(0.0) 


The  origin  Z0  in  the  polynomial 
expansion  for  the  vector  potential. 
A  (xty)=re(sum  Cn* ( Z-Z0 ) **n ) ,  the 
derivatives  of  which  give  the  field 
and  gradient.  Note  -  for 
cylindrical  symmetry  problems, 
XORG=0,  always. 
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CON(42)=KMIN 

(1) 

CON(43)*KTOP 

(UMAX) 

CON (44 ) =LMIN 

(1) 

C0N(45)=LT0P 

(1) 


CON(46)=ITYPB 

(2) 


CON ( 47 ) =  W2ND 
(0.125) 


The  (K,L)  liaits  of  the  region  in 
which  the  fields  and  gradients  are 
to  be  calculated  at  each  aesh 
point.  See  CON(54)  through  CON(57) 
for  an  alternative  scheae  of 
defining  the  region  in  which  the 
fields  and  gradients  are  to  be 
calculated. 


A  code  specifying  the  problea 
syaaetry  and  deteraining  which  "C  ” 
tera  appears  in  F(s),  and  whether 
the  HC  H  are  real,  iaaginary,  or 
coaplex ? 

ITYPB=1,  No  syaaetry. 

1TYPE=2,  Midplane  syaaetry 
ITYPE=3,  Elliptical  quad  aperture 
ITYPB=4,  Syaaetrical  quad. 

ITYPB=5,  Skew  elliptical  quad 
aperture. 

ITYPB=6,  Syaaetrical  "H"  aagnet  or 
elliptical  aextupole  aperture. 
ITYPB=7,  Syaaetrical  sextupole 
ITYPE=8,  Elliptical  octupole 

aperture. 

ITYPB=9,  Syaaetrical  octupole. 

For  all  of  the  above  syaaetry 
codes,  except  for  ITYPB=1,  or  =5, 
field  lines  are  perpendicular  to 
the  x-axis.  For  ITYPB=5,  the  x- 
axis  is  a  field  line. 

For  cylindrical  syaaetry  (vector) 
prob leas : 

I TYPB  =  1 ,  No  syaaetry 

ITYPE=2,  Midplane  syaaetry,  r-axis 

has  field  lines  perpendicular. 

For  cylindrical  syaaetry  (scalar) 
probleas , 

ITYPE=1,  No  syaaetry 

ITYPE=2,  Midplane  syaaetry,  r-axis 

has  V  lines  perpendicular. 

ITYPB=3,  Midplane  syaaetry,  r-axis 
is  a  V=constant  line. 


The  weight  factor  for  the  second 
nearest  neighbors  used  in 
deteraining  the  "C  "  in  the 
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expansion  for  A(x,y). 


CON(48) = ISECND  Indicator  to  use  second  neerest 

neighbors  in  determining  the  "C  ” 
in  the  expansion  for  A(x,y). 

ISECND= 1 ,  Use  first  and  second 
neighbors 

ISECND=0,  Use  first  neighbors  only. 


CON ( 54 ) =XMIN 

(0.0) 

CON (55 ) =XMAX 

(0.0) 

CON ( 56) =  YMIN 

(0.0) 

CON ( 57 ) = YMAX 

(0.0) 


The  (x,y)  limits  of  a  grid  over 
which  the  fields  and  gradients  are 
to  be  calculated.  Delta-x  and 
Delta-y  in  this  grid  are  specified 
by  CON ( 43 ) =  KT0P  and  CON(45)=LTOP 
as  follows: 

Dx= ( XMAX-XMIN ) / ( KTOP- 1 )  and 
Dy= (YMAX- YMIN) /(LTOP-1) . 


4.2. le  Current  Parameters  and  Options 


CON ( 8 ) =BDES 
(1.0E15) 


BDES= : B :  at  the  location 
KBZ8R0=C0N ( 40 ) ,  LBZERO=CON ( 4 1 ) .  If 
BDES  is  not  equal  to  1.0E15,  the 
current  factor,  XJFACT=CON( 66 ) , 
will  be  adjusted  so  that  :BI=BDES 
within  a  tolerance  XJTOL=CON ( 67 ) . 


CON ( 40 )  =  KBZERO 

(1) 

CON (41 )=LBZERO 

(1) 


The  (K,L)  coordinates  specifying 
the  location  of  BDES  for  adjusting 
the  current  factor. 


CON(49) =NFIL 

(0) 


The  number  of  current  filaments  to 
be  read  in  as  data  by  problem 
POISSON. 


CON (66)= XJ FACT 

(1.0) 


The  factor  by  which 
densities  (except 

filaments)  will  be 
program  POISSON, 
indicates  a  scalar 


all  current 
current 
scaled  in 
XJFACT=0. 0, 
potential 


problem  (no  currents).  LATTICE 
Mill  automatically  set  XJFACT=0.0 
when  IB0UND=-1 . 


4.2. If  Over  Relaxation  Factors 


CON( 74 ) =RH0PT1 
(1.90) 


See  RHOAIR=CON(75) 


CON(75)=RHOAIR 

(1.90) 


The  over-relaxation  factor  in 
POISSON  for  air  and  interface 
points  (and  iron  points  during  a 
au-f inite-but  constant  solution). 
This  factor  ia  automatically 
optimized  during  the  iteration  if 
the  initial  value=RHOPTl=CON(74) . 


CON (77 ) =RHOFB 

(1.0) 


The  over-relaxation  factor  in 
POISSON  for  iron  points  during  a 
au-f inite-variable  solution. 


CON ( 78) -RHOGAM 
(0.08) 


The  under-relaxation  factor  in 
POISSON  for  gamma  during  a  au- 
f inite-variable  solution. 


CON(80) = ISKIP 

(1) 


The  number  of 
recalculating  the 
au-f inite-variable 


cycles  between 
gammas  during  a 
solution. 


$.?■ ll _ Convergence  Criteria 


CON (85)=BPSILA 
(5. OB-7) 


The  convergence  criterion  for  the 
potential  solution  of  air  and 
interface  points  (and  for  iron 
points  during  a  au-f inite-but 
constant  solution). 
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Sec.  4.2.2  Paraaeters  for  Space-Charge  Flow 


This  section  describes  the  CON  s  which  are  used  in  POISSON  for 
the  space-charge  flow  problem.  The  description  of  the  CON’s  for 
space-charge  flow  are  described  in  Sec.  3.7.2. 


Sec.  4.2.2a  Printing  Options  for  the  Results  of  POISSON 


MPRINT  (CON(143)),  controls  the  printing  options  for  the  results 
of  POISSON.  There  are  two  output  files  where  the  results  are 


printed: 


1)  POISSON. OUT 

2)  POSITION. OUT 


1)  "POISSON.  OUT"  contains  the  printed  output  of  the  CON’s, 
cathode  data,  potential  at  each  mesh  point,  the  charge  assigned 
to  each  mesh  point,  current  emitted,  convergence  information, 
perveance  and  total  beam  current. 


MPRINT  has  three  levels  that  may  be  specified  in  order  to 
retrieve  some  or  all  of  the  above  information.  These  three 
levels  are: 


MPRINT= 1  When  MPRINT  is  set  to  one,  the  following  items 
are  printed: 

1)  CON’s 

2)  Number  of  iterations  for  the  potential  solution  to 
converge  for  each  current  cycle. 

-  77  - 


few* 


•  ,> 
a  '  »  >  « 

m 

» 

li 


%■  •.*  >, 


3)  Total  beam  current,  perveance  and  DCBAR  for  each 
current  cycle. 


m 


MPRI WT=2  In  addition  to  the  above  iteas,  the  following 


iteas  are  printed  when  MPRINT  is  set  to  two: 


1)  cathode  data 


2)  CBAR  and  eaitted  current  of  each  eaitting  point  for 


each  current  cycle. 


MPRINT=3  In  addition  to  the  above  iteas,  the  following 


iteas  are  printed  when  MPRINT  is  set  to  three: 


1)  the  potential  at  each  aesh  point  for  each  current 


cycle. 


2)  charge  assigned  to  each  aesh  point  for  each  current 


cycle. 


2)  "POSITION. OUT”  contains  the  trajectory  inforaation  of  the 


electrons  when  MPRINT  is  greater  than  one.  Tiae  increaent,  z(t). 


r(t),  z-velocity,  r-velocity,  theta-velocity,  perpendicular 


velocity  and  two  relativistic  gaaaas  are  printed  at 


INCPLT=CON( 146)  increments.  For  problems  which  have  a  large 


aaount  of  tiae  steps,  INCPLT  should  be  set  accordingly  so  the 


plotting  arrays  do  not  overflow.  A  value  of  INCPLT=3  is  typical 


for  problems  with  many  tiae  steps. 


Two  gaaaas  are  calculated  by  Eq.  (5.3)  and  Bq.  (5.25)  of  the 


thesis.  These  two  gaaaas  should  alaost  equal  each  other.  If 


they  do  not,  energy  is  not  being  conserved,  and  the  tiae  step  is 


probably  too  large. 
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Sec.  4.3  Computation  of  the  Electron  Trajectories  of  MODGUN 


In  Sec.  3.8,  LATTICE  created  the  triangular  nesh  and  created  the 
data  file  POISSON.DAT  (  duap  0  )  for  the  problem  MODGUN.  POISSON 
will  read  in  POISSON.DAT  and  compute  the  electron  trajectories. 


Once  POISSON  has  been  executed,  the  program  will  prompt  on  the 
screen : 


INPUT  DUMP  NUMBER,  "NUM* 


NUM  is  the  "dump"  number  to  be  read  in,  as  discussed  in  Sec.  4.1. 
Dump  0  is  to  be  read  in  by  POISSON,  therefore  the  user  responds 
by  typing  in  0  and  then  a  return. 


POISSON  will  next  ask  the  user  to: 


INPUT  THE  CONS  YOU  WANT  TO  CHANGE 


As  discussed  in  Sec.  4.1,  this  is  the  second  group  of  input  data. 
At  this  point  the  user  may  want  to  change  some  CON's  that  were 
not  specified  in  LATTICE.  For  example,  to  get  a  full  diagnostic 
print  out,  MPRINT  is  set  to  three.  The  following  is  typed  on  th« 
terminal : 


*143  3  S 


where  the  format  of  this  line  is  discussed  in  Sec.  3.1 
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The  prograa  PLOT,  plots  the  aesh,  equipotent ial  lines  and 


electron  trajectories.  PLOT  uses  Hewlett  Packard's  Advanced 
Graphics  Package  (AGP)  to  perfora  the  plotting. 


A  ainiaua  of  1  input  line  aust  be  provided  for  PLOT: 

1)  The  first  input  line  gives  values  for  five  paraaeters  for 
the  plot: 


"NUM"  is  the  "duap"  number  to  be  read.  The  aesh  is  written 
on  duaps  "0”  and  "1",  and  the  equipotent ial  lines  and 
electron  trajectories  are  written  on  dump  "1”. 

" ITRI"  is  a  flag  whether  or  not  to  plot  the  triangles  of  the 
aesh. 

ITRI=0,  do  not  plot  the  triangles 
ITRI=1,  plot  the  triangles 


"NPHI "  is  the  nuaber  of  equipotent i al 
The  prograa  does  not  plot  the  ainiaua 
values  which  is  usually  a  boundary. 


lines  to  be 
or  maximum 


plotted . 
potential 


"I NAP"  is  the 
ainiaua  and 
equipotent ial 


indicator  to  read,  on  the  next  data  card, 
aaxinum  values  (AMIN  and  AMAX)  of 
lines  to  be  plotted. 


the 

the 
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INAP=0,  do  not  read  AMIN  and  AMAX 

INAP= 1 ,  read  AMIN  and  AMAX  on  the  next  input  line 


The  valuea  plotted  are  between  (AMAX-DELTA)  and  ( AMIN+DELTA) 
where 


DELTA= ( AMAX-AMIN) / (NPHI+1 ) 

" IRAYS"  ia  a  flag  whether  or  not  to  plot  the  electron 
trajectoriea. 

IRAYS=0,  do  not  plot  the  trajectoriea 
IRAYS=1,  plot  the  trajectoriea 

2)  The  aecond  input  line  is  optional  and  gives  values  for  four 
paraaeters  for  the  plot:  XMIN,  XMAX,  YMIN ,  YMAX. 

These  four  paraaeters  specify  the  liaits  for  the  plot,  which 
aay  be  any  part  of  the  problea.  If  this  data  is  skipped, 
the  entire  problea  will  be  plotted. 

Other  options  such  as  sending  the  plot  to  the  printer  and  zooaing 
in  on  the  plot  are  also  available  by  answering  the  appropriate 
questions  proapted  by  PLOT. 

Sec.  6.2  Plotting  the  Results  of  MODGUN 

The  triangular  aesh  aay  be  plotted  after  LATTICE  has  been 
executed.  The  aesh  of  MODGUN,  Fig.  3.12,  was  plotted  by 
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